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Abstract 



We point out that the measurement of the L/E distribution of atmospheric neutrino events 
with sufficient resolution in order to observe the characteristic modulations of neutrino oscillations 
would provide a unique way to test the — ► v x oscillations nature of the atmospheric neutrino 
anomaly jl) in the parameter range 5 x 10~ 4 < Am 2 < 2 x 1CP 3 eV 2 . With an integrated exposure of 
a few tens of kiloton-year, this method can cover the lower part of the Am 2 region of the currently 
ON ■ allowed atmospheric neutrino solution and is therefore complementary to the one explorable with 

planned long-baseline neutrino beam experiments (e.g. Refs. || f|). We stress that experimentally 
this method relies on the ability to measure high energy muons (E^ > 5 GeV) . A large magnetized 
calorimeter of several kilotons would be ideally suited for this measurement. The existing LEP 
detectors augmented with atmospheric neutrino triggers could also be used to perform such a 
<3\ ' study. 

^ \ 

^ 1 Introduction 

Recent SuperKamiokande data show an asymmetry in the rate of upward and downward atmo- 
spheric muon-like events (l], In order to verify that this phenomenon is really due to neutrino 
oscillations, an effective method consists in observing the sin 2 (1.27Am 2 L/i?) modulation charac- 
teristic of a given Am 2 when the event rate is plotted as a function of the reconstructed L/E of the 
events. Such a method, capable of measuring Am 2 in the range between 5 x 10" 4 and 2 x 10" 3 eV 2 
exploiting atmospheric neutrino events, has been advocated in a recent note||. The method has 
the advantage of being practically insensitive to the precise knowledge of the atmospheric neutrino 
flux since the oscillation pattern is found by striking dips in the L/E distribution while the neutrino 
interaction spectrum is known to be a slowly varying function of L/E. It was shown that the main 
experimental requirement is that the error on L/E is the following: 

5{L/E) < tt/(1.27Ato 2 ) (1) 

This implies that precise measurements of energy and direction of both the muon and hadrons 
are needed in order to reconstruct precisely the neutrino L/E. Note that all past and present 
large scale underground atmospheric neutrino experiments are unable to determine accurately the 
momentum of high energy muons (e.g. E^ > 5 GeV). This implies that the capability to study 
the L/E distribution is currently quite limited. The ability to measure high momentum muons 
requires the construction of a magnetized calorimeter with a mass of several kilotons in order to 
have enough of an event interaction rate. 

In section |L we define a set of selection cuts that improve the reconstruction of the L/E 
distribution of events which do not strongly rely on the hadronic energy resolution. In section [|, 
we show how the oscillation can be found and the oscillation parameters determined in a generic 
magnetized calorimeter. In section^, we argue that the existing LEP detectors — able to measure 
muon momentum and direction with a coarse hadronic energy resolution could fulfill the required 
criteria to observe the L/E modulations. 
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2 Optimizing the L/E resolution 



The atmospheric neutrino interaction energy spectrum for E > 1 GeV has roughly a E 
dependence ||. The basic idea to satisfy Eq. [I] consists in selecting only the neutrino events where 
the muon direction is close to that of the incoming neutrino. Given the energy dependence of the 
atmospheric neutrinos, the events with high energy muons are more likely to fulfill the condition. 
Therefore, there is a cut on the muon momentum for which the total event rate and the L/E 
resolution are optimized. Moreover the effective cut in the neutrino energy removes the smearing 
due to the Fermi motion of the target nucleon. 

On a sample of simulated atmospheric neutrino interactions, we placed a cut on the muon energy 
at 2.5 GeV to select a sample of events where an 80% fraction has the opening angle smaller that 
10°. The remaining 20% is mainly due to high-y high neutrino energy events and is small due to the 
E~ 1,7 dependence of the neutrino interaction spectrum. It can be further reduced requiring small 
hadronic energy deposition; a cut of Ei ia d < 0.5 x satisfies the above requirements. We point 
out that for this purpose a good hadronic energy resolution is not needed. The neutrino energy 
is reconstructed using both muon and hadron energy but the neutrino energy resolution is only 
partially affected by the rough hadronic energy resolution because of the applied selection criteria. 

Another important cut is applied on the incoming neutrino zenith angle 9\ in order to minimize 
the error on L 1 only the events with 9 < 75° or 9 > 105° should be retained (see Ref. to well 
define the sets of upward and downward events. The selection efficiency of this method is close to 
15 % for events in the detector with neutrino energy above 1.0 GeV. The selection efficiencies for 
neutrinos and antineutrinos are summarized in Table |[ The efficiency for antineutrinos is higher 
due to da/dy oc (1 — y) 2 compared to the flat da/dy for neutrinos. 



3 Oscillation parameters determination 

To illustrate the method, we consider a generic calorimetric detector with muon identification and 
measurement capability. The assumed resolutions are listed in Table |^. 

Given the oscillation parameter Am 2 of interest, we expect that the downward atmospheric 
neutrinos above 2.5 GeV will not have oscillated while the upward going neutrinos will show the 
oscillation^). The sensitivity to oscillations is therefore obtained by comparing the measured L/E 
distributions of upward and downward going neutrino events. For upward going events, the distance 

L is obtained from the cos 9 with the expression L = ^R 2 - (R - d) 2 sin 2 9 -(R-d) cos 9 where 
R is approximately the earth radius and d is the depth of the atmosphere. For downward going 
events, the mirror- distance L' is calculated replacing in the expression 9 — ► —6, in order obtain an 
"unoscillated" upward going V /E spectrum. This method only relies on the fact that the angular 
distribution of atmospheric neutrinos is symmetric in cos 9 for energies above 5 GeV and below 
5 GeV the geomagnetic effect and the solar activity perturbances are small Q. In the considered 
cos 9 region (see section ^|), the geomagnetic correction is even smaller. The — ► v x conversions 
are extracted by searching for the characteristic pattern of oscillations in the L/E distribution 
normalized by the L'/E distribution (R = (L/E)/(L'/E)). 

In Figures |^J^ we show the measured L/E distribution in presence of — > v x oscillations with 
parameters Am 2 = 1 x 10~ 3 and Am 2 — 5 x 10~ 4 and sin 2 29 = 0.9 after an exposure of 30 
kt-year for upward going and downward going muons. The oscillation parameters are extracted 
by comparing the two distributions. For simplicity, we show the ratio R of the two distributions 



Selection 


Fraction(%) 


Fraction(%) 


S(L/E) 
RMS (Km/GcV) 


All {E v > 1 GeV) 


100 


100 


765 


E^ > 2.5 GeV 


31.7 


45.1 


560 


E had < 0.5£ M 


13.7 


27.0 


350 


9 < 75° or 9 > 105° 


11.2 


21.9 


345 



Table 1: Neutrino and anti- neutrino event selections. 
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Muons 



Hadrons 



Generic calorimeter: 



SE/E 
56 



5% 80% /y/E 



10 mrad not used 



Table 2: Resolutions assumed for a generic calorimetric detector with muon identification and mea- 
surement capability. 

where we have assumed symmetric errors. The position of the first minimum in the R distribution 
(corresponding to an (L / ' E) m i n ) yields an estimate of the oscillation parameter Am 2 through the 
oscillation formula. The best estimate of the parameters in the (Am 2 , sin 2 26) plane are obtained 
by fitting the oscillation formula to the observed R distribution. The confidence levels (see Figures 
for the determination of the mass difference and mixing angle at the 90, 95 and 99% shows 
the power of this method. 

4 Atmospheric Neutrinos in LEP detectors 

The existing LEP detectors ( ALEPH [§, DELPHI g, L3@ and OPAL0), equipped with some 
additional apparata, could fulfill the required criteria; namely: 

• large mass (several kt); 

• magnetized muon spectrometer; 

• segmented hadron calorimeter; 

• very precise tracking chambers (TPC's and L3 MUCH). 

We note that for our purpose the three detectors ALEPH, DELPHI and OPAL have very similar 
characteristics with most of their mass concentrated in their hadronic calorimeters acting as return 
yokes for their inner magnetic field. The L3 detector has a large mass in its iron magnet surrounding 
completely the large precise muon chambers and the inner sub-detectors. We estimate that a useful 
mass of « 3 kt can be obtained with the iron of the L3 detector magnet; the hadron calorimeters 
of ALEPH, OPAL and DELPHI amount to 9 kt; about one more kt can be added if one includes 
the electromagnetic calorimeters. 

For the events occurring in the hadronic calorimeters the neutrino energy is estimated as the 
total visible energy and an appropriate cut on the hadronic activity will be used to improve the 
correlation between neutrino and muon direction. We stress the importance of the magnetic field 
available in all LEP detectors in order to obtain a reliable muon momentum measurement: 

• In the case of ALEPH, the acceptance for the muon to traverse the inner tracking detector 
is large since the hadron calorimeter was built to completely cover the inner tracker. In this 
case a very precise determination of the muon momentum is available; for the events with 
this topology the measurement error is dominated by the error of the backward extrapolation 
into the hadron calorimeter. For the muons that do not cross the inner tracker, the fine 
segmentation of the hadron calorimeter should permit the momentum to be measured by 
curvature (using the residual magnetic field) or by range-out. 

• For the events in the L3 magnet, the hadronic energy cannot be reliably estimated. Never- 
theless, we recall that the cut on muon momentum suppresses topologies with large hadronic 
energy and a further reduction can be obtained by requiring no other activity that the muon 
in the muon chambers. So we expect that the muon momentum vector will be a good es- 
timator of the incoming neutrino energy and direction. To efficiently veto cosmic ray muon 
backgrounds, we expect that the L3 magnet will have to be covered with active large area 
detectors (e.g. RPC's). 

We deduce from the technical reports of the LEP experiments ^, ||, || [hJ that the resolutions 
for reconstructing the events in the topologies described above are well matched with the one 
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Figure 1: Measured L/E distribution in presence of — > v x oscillations with parameters Am 2 = 
1 x 10 -3 and sin 2 20 = 0.9 after an exposure of 30 kt-year for a) upward going muons; b) downward 
going muons; c) the ratio R; d) the confidence levels for the determination of the mass and mixing 
angle. 



hypothesized for the generic calorimeter of section A more precise description and discussion on 
the use of the LEP detectors will be given after detailed simulations will have been performed. 



5 Conclusions 

In this note we have shown that the neutrino oscillation parameter, Am 2 , can be precisely deter- 
mined through the measure of atmospheric neutrino flux as a function of L/E in the search region 
between 5 x 1CP 4 and 2 x 10 -3 eV 2 even with a large mass detector with rough hadronic resolu- 
tion in the GeV range. The four LEP detectors have a total mass of about 10 kt and resolutions 
matching the requirements to use the L/E method. It could be very interesting to explore their 
capabilities as atmospheric neutrino detectors. 
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